Non-Fermi liquid fixed point for an imbalanced gas of fermions in 1 + e dimensions 
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We consider a gas of two species of fermions with population imbalance. Using the renormalization group in 
d = 1 + e dimensions, we show that for spinless fermions and e > a fixed point appears at finite attractive 
coupling where the quasiparticle residue vanishes, and identify this with the transition to Larkin-Ovchinnikov- 
Fulde-Ferrell order (inhomogeneous superconductivity). When the two species of fermions also carry spin 
degrees of freedom we find a fixed point indicating a transition to spin density wave order. 

PACS numbers: 



Experiments on ultracold atomic gases allow fermionic 
pairing phenomena to be investigated with a degree of con- 
trol and purity hitherto unknown in solid state systems. 
The preeminent example is the observation of the crossover 
from Bose-Einstein condensation (BEC) to Bardeen-Cooper- 
Schreiffer (BCS) superfluidity effected by tuning the scatter- 
ing length between two atomic components using a Feshbach 
resonance fl] [3] g) . 

It is standard lore that pairing between two components at 
equal densities - hence with equal Fermi wavevectors - occurs 
for arbitrarily weak interactions in the ground state, though the 
transition temperature may become very small. By contrast an 
imbalance in the two populations requires a sufficiently strong 
attraction before pairing takes place [5 1. The first experiments 
on the imbalanced system ||6] [7] [8] |9) confirmed this picture, 
together with the expectation that the transition between the 
normal and BEC superfluid states is first order, leading to 
phase separation. 

An alternative route to pairing in the imbalanced case 
was introduced in Refs. ifTUl ITT1 . The Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) state, as it is now called, is formed from 
pairs with center-of-mass momentum equal to the difference 
in Fermi momenta. Within mean-field theory the LOFF phase 
occupies a rather small part of the phase diagram in terms of 
imbalance and interaction strength [12] (see Fig. [T]>. In lower 
spatial dimension, the LOFF phase is expected to be more 
prominent [13]. Experimentally, there is some evidence that 
the LOFF state occurs in the heavy fermion superconductor 
CeCoIn 5 lfl4l[T5l. 

In distinction to the normal-BCS transition, the normal- 
LOFF transition is expected to be continuous at low temper- 
atures. It is therefore somewhat surprising that to date there 
has been no attempt to understand the nature of the quantum 
phase transition out of the Fermi liquid state that occurs with 
increasingly attractive interaction. It is the purpose of this Let- 
ter to provide that understanding. 

Our main findings can be summarized as follows. Along 
the line that forms the phase boundary between the normal 
and LOFF states in a diagram of polarization versus inter- 
action strength (solid red line in Fig. [TJ the critical state is 
characterized by a singular interaction between fermions of 
the two species with anti-parallel momenta. This leads to an 




FIG. 1 : (Colour online) Schematic phase diagram showing the con- 
tinuous quantum phase transition between the normal Fermi liquid 
and LOFF states as a function of interaction strength and imbalance. 
The LOFF-BCS transition is first order, leading to phase separation. 
Inset: vanishing of the quasiparticle residue in the Fermi liquid as the 
transition is approached along the bold horizontal arrow. Note that 
the exponent varies continuously with V '. 



incoherent spectral function for particles at the Fermi surface 

A S (\K\ = K FtS ,w) (1) 

where s = a,b labels the species, with an exponent r\ that 
varies continuously along the critical line as a function of the 
polarization V = (n a — rib)/(n a + n&). This line is a higher 
dimensional analog of the Luttinger liquid in one dimension. 
For weak polarization, V <C 1, this result is valid for spatial 
dimension 1 < d < 2. As the transition is approached from 
the normal phase, the quasiparticle residue vanishes contin- 
uously. It is our hope that some subset of these predictions 
can be probed in an ultracold gas by the recently developed 
technique of momentum-resolved RF spectroscopy |[T6l [T71 . 
Towards the end of this work we present an extension to the 
problem of spin density wave ordering in systems with par- 
ticle and hole Fermi surfaces, which may of relevance to the 
recently discovered iron-based pnictide superconductors ifTHl . 

In what follows we assume n a > rib an d hence the Fermi 
wave vectors obey Kp, a > Kp,b- The energy of a fermion of 
species s is given by £ S (K) = e s (K) — fi s , where e s (K) = 
K 2 /2m s and /i s = Kp s /2m s , with mass m s . Near their 
respective Fermi surfaces, where we expect the important 
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components and where we have set h = 1. We start by moti- 
vating the problem of pairing with unequal Fermi surfaces. 

Cooper's problem Consider a pair of fermions of either 
species, above their respective Fermi seas. The two fermions 
interact only with each other, with the Fermi seas serving only 
to block states below the Fermi level [ 19 1. In vacuum the two- 
particle Schrodinger equation takes the form 



2m a 2m b 



V ' ' V5(r a -T b )-E 



4>(r a ,r b ) = (3) 



After imposing the restrictions due to Pauli blocking, Eq. ([3]) 
is equivalent to the condition 



FIG. 2: (Colour online) Pairing of majority (blue) and minority (red) 
fermions at fixed total momentum Q. 



physics to occur, the fermions have an approximately linear 
dispersion: £ S (K) = vp, s k + 0(k 2 ), where vp tS = Kp. s /m s 
is the Fermi velocity and k = K — Kp. s is the momentum 
relative to the surface. The effective fermionic Hamiltonian 
we will work with is 



W = £k(K)$(K)^(K) 



K,s 



+ V ]T ^(Q-K)^(K)^ a (K0V6(Q-K'), (2) 

KK'Q 

with a point interaction V that acts only between the different 
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where Q is the center-of-mass momentum. A bound state 
corresponds to E = n a + n b +E b for some E b < 0, and to find 
such a solution for V small requires that the momenta of the 
two particles are close to antiparallel. Setting |Q| = Kp^ a — 
Kpj, and imposing a momentum shell cut-off of thickness A 
around each Fermi surface, we see that the angle between p 
and Q is limited as p — ATp,b — > by the condition (see Fig. [2]) 
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Setting E b — to find the threshold for bound state formation, 
the right hand side of Eq. Q becomes for A <C Kp s 
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where S d = 2ir d / 2 /T(d/2) is the area of the unit sphere in d 
dimensions and vp — (vp^ a + vp^ b )/2. For small e = d — 1, 
S e ~ e, Eq. Q becomes 

1 (2Kp. a K F . b y /2 K^-^ 

, (6) 
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where we have introduced a small infrared cutoff /i. In the 
limit e — > there is a logarithmic singularity log(A//i), lead- 
ing to a logarithmically small bound state energy for arbi- 
trary negative V. This divergence in the particle-particle (PP) 
scattering channel is a consequence of Fermi surface nesting, 
£a(Q ~ P) + £f>(p) = 0, for antiparallel fermions of the two 
species, familiar from the usual Cooper problem and the BCS 
theory. The feature that we wish to emphasize is that in the 
presence of imbalanced Fermi surfaces, the interaction of an- 
tiparallel fermions only produces a logarithm in d = 1. In 
the renormalization group (RG) sense it is marginal in d = 1, 



becoming irrelevant for d > 1. The above calculation ne- 
glects however the particle-hole (PH) contribution to scatter- 
ing, which motivates the following more careful RG analysis. 

RG Calculation The scattering behaviour is encapsulated 
by the four point vertex function, T, the only contributions 
to which at one loop order are the bubble diagrams shown in 
Fig. [3] corresponding to PP and PH excitations. For pair mo- 
mentum Q = Kp a — Kp b and vanishing external frequency, 
uj — > 0, these may be evaluated for small angles. The result 
for the PP bubble is 
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Here the factor of two relative to Eq. (|6]l is due to an equal 
contribution from pairing below the Fermi surfaces. A compa- 
rable treatment results in a similar expression for the PH bub- 
ble, but with the opposite sign and Q — > Q' = Kp a + Kp tb . 
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FIG. 3: The two diagrams contributing to the vertex at one loop order. 
Dashed and solid lines indicate impurity (a) and gas (b) propagators 
respectively. The notation is explained in the text. 
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FIG. 4: Flow diagram for the spinless case, indicating Fermi liquid 
(FL) and LOFF phases. 



Combining the terms, we have for scattering of an antiparallel 
pair at the Fermi surface 



r = v -v 



2 2S e (2^K Ffi K F , a y' 2 
(2 7 r) 1 +^ F e 2 



A £ / 2 F 



with 



F = 



Q- £ / 2 - Q'- £ / 2 



(8) 



(9) 



We see that there is cancellation of the two channels if d = 1 
or Kpjj/Kp.a = 0, so that F = 0. 

By demanding cut-off independence, A^iT = and defin- 
ing a dimensionless coupling, g, via 



S e (2Ay/K F , b K F , a ) 
(2ir) 1+e ev F 



e/2 



V, 



we obtain the beta function 



(10) 



(11) 



Hence there is a non-trivial fixed point = —jp, which is 
unstable as illustrated in Fig. |4] Notice the difference from the 
more familiar situation, typified by the Wilson-Fisher fixed 
point, in which the interacting and free fixed points merge as 
the critical dimension is approached. In this case, the vanish- 
ing of F as e — > means that (3 vanishes, leaving the fixed 
point at finite coupling. The e = cancellation extends to all 
orders (this may be shown, for example, using Ward identi- 
ties EOl ) and hence (3(g) = for a range of g, indicating the 
existence of the well known Luttinger liquid critical phase. 

For e > and weak imbalance, Kp ia ~ K F:b , the PP con- 
tribution to Eq. ([9]) dominates and the fixed point occurs at 



weak coupling 
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Integrating the flow for \g\ <C |<?*| yields g = go(A/ Ao) e / 2 , 
where go and Ao are the initial conditions for the scaling. Con- 
sequently V is independent of A. This is the expected result 
for couplings when V = and Q ^ lETl . In contrast, when 
g = we find V ~ — A~ e / 2 so that the coupling becomes in- 
creasingly attractive as A decreases. The situation for V -C 1 
corresponds to a microscopic perturbation theory in which the 
vertex is given as a geometric sum of PP bubbles and the re- 
sulting expression can be extended to e = 1. With this in mind 
we expect the weak polarization fixed point result to be valid 
in d = 2. 

For strong imbalance we may expand F as a power series 
in Kp.bl K F:a , F w eK F ^/ K Fya . As a result j3 is given by 
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The fixed point is now at strong coupling, independent of e, 
and higher order contributions should then be taken into ac- 
count. As such it is not possible to make quantitative state- 
ments about the strongly imbalanced case, V ~ 1 using this 
method. 

We now turn to the fermion self energy. The one loop 
contribution produces a routine frequency independent shift 
which may be absorbed by the chemical potential, so as to 
keep K F s fixed. The first significant behaviour is found by 
calculating the two loop 'rising sun' diagram, yielding 



0ReE?(|K| =K Fi s,lj) 



9 2 C s , t 
4e 



0(g 3 ). (14) 



Though the coefficient C Si£ depends on the ratio K F ^ a / K F ^, 
it goes to unity in the limit e — * 0. Note that Eq. ( p"4| ) 
arises from scattering of almost antiparallel particles of the 
two species. 

A independence of the physical correlation functions ne- 
cessitates the introduction of a multiplicative field renormal- 
ization ^P hys = Z 1 / 2 ip s . Z is fixed by demanding A indepen- 
dence of the physical Greens function 



G^ S (\K\^K F ^) 



-£*(|K|=A>,„, W )- 



This implies 



d\nZ q „, . q „, 

d^A = -i^ = h + °^ 



(15) 



(16) 



so that at g = g±, Z — * as A — * 0. At the fixed point the 
scaling behaviour of the propagator is given by 



(17) 



where the critical exponent r\ — g 2 /^ = e 2 /16F 2 . Hence the 
quasiparticle pole is replaced by a branch cut with attendant 
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power law divergent spectral function. It is also enlightening 
to determine the dependence of Z on the coupling as it ap- 
proaches g+. By linearizing the flow equations near the fixed 
point we may derive the following relation between the start 
and end points of the flow, (Zq, go) and (Z,g) respectively, 



Z= Zo 



9o - 9 



e/8F 2 



exp 



F 



(g - go) 



(18) 



Using the above we may extract the value of the Fermi liquid 
quasiparticle residue, Z, due to an initial coupling go. If we 
take go > g* so that the coupling flows to g = 0, the final 
value of Z will depend on go as Z ~ (go — g+) e > 8F (see inset 
to Fig. |T}. 

It is possible to extend the calculation above to treat com- 
ponents that themselves have internal degrees of freedom. In 
the case that a and b are both spin half particles, the quartic 
part of the Hamiltonian becomes 

'Hint = [(V e ||<J ff)<7 > + Ve^<j-a')^\^\^,^a,o^b,o' 

(7,(7' 

(19) 

where we have supressed momentum labels. Clearly there is 
no discernible difference between processes due to V e \\ and 
those due to Vrfii . Therefore we define Vi| = V^n +Vd|| - If Vri = 



Vd± Eq. ( 19 1 is equivalent to an anisotropic Heisenberg in- 



teraction, with J z — 4 Vfi and J xy = 2V e ±. Keeping only the 
lowest order terms (so that F = 0) and defining dimension- 
less couplings, gi — ViAe, where A — (Kp^Ay/ 2 /nvpe, 
one finds 



P(g\\) 

0{9dx) 
(3{9ex) 



:9d± 



- 2 



~2 9e± ~ 9e±(9d± ~ 9\\)- 



(20) 
(21) 
(22) 



The first two equations are equivalent if Vfi = — Vd_i_. Focus- 
ing on this case we discern that there is a non-trivial fixed 
point at (g\u, g e ±*) = (e/4, e/2). Concomitantly, and in 
terms of the Heisenberg exchange couplings, there is a transi- 
tion to spin density wave order at antiferromagnetic exchange 



Jz* = Jxy* = envp/ (Kp jb A) 



;/2 



An experimentally relevant scenario is that of interacting 
electron and hole pockets [22 1. Such systems have been inves- 
tigated extensively in the context of excitonic transitions and 
antiferromagnetism (see for example E3l ) and recently be- 
cause of interest in the high temperature superconductivity of 
iron pnictides ll24l . In the calculation presented here, switch- 
ing one of the components from particle-like to hole-like al- 
ters the dispersion as £k.s — * £k,/i = —vp y hk and changes 
the pairing wavevector to |Q| = Kp e + Kph where the 
subscripts e and h indicate electrons and holes respectively. 



In terms of Eqs. (20 - 22 1 the effect is to flip the sign of the 
0(g 2 ) terms. Again taking Vfi = — Vd±, the fixed point is at 
(g\\*,g e ±*) = (— e/4, e/2) or equivalently, for ferromagnetic 
exchange J zi , = -J xy * = -eitvp / (Kp^Af/ 2 . 

In this letter we have applied the renormalization group to 
polarized two species fermi gases in 1 + e dimensions. The 
central result for spinless fermions is the appearance of a non- 
Fermi liquid fixed point, characterised by finite pairing wave 
vector (centre of mass momentum) and power law spectral 
function, with an exponent that depends on the polarization. 
For fermions that carry spin we have shown that there are non- 
trivial fixed points that describe magnetic ordering of the spin 
density wave type. AL gratefully acknowledges the support 
of the NSF under grant DMR-0846788. 
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